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We compute the average shape of trajectories of some one-dimensional stochastic processes x{t) 
in the {t, x) plane during an excursion, i.e. between two successive returns to a reference value, find- 
• ing that it obeys a scaling form. For uncorrelated random walks the average shape is semicircular, 

independently from the single increments distribution, as long as it is symmetric. Such universality 
Mh , extends to biased random walks and Levy flights, with the exception of a particular class of biased 

Levy flights. Adding a linear damping term destroys scaling and leads asymptotically to flat excur- 
sions. The introduction of short and long ranged noise correlations induces non trivial asymmetric 
shapes, which are studied numerically. 

rS^ • PACS numbers: 05.40.-a,75.60.Ej,05.45.Tp 
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^ ; I. INTRODUCTION 

> ' 

Many disordered systems respond to external solicitations by producing noise with power law features, that can be 
^2 ' modeled in terms of avalanches. A notable example of such phenomena is the Barkhausen effect, first observed about 
a century ago by recording the noise produced by the reversal of large domains in a ferromagnet. The Barkhausen 
' ' noise has been incessantly investigated, both because of its practical application as a nondestructive method to test 
magnetic materials, and because of its conceptual relevance for the understanding of the magnetization dynamics 
on a microscopic scale Experiments show that both the size and the duration of avalanches of spin reversal are 
^ power law distributed over several decades. The exponents characterizing these power laws are often used to identify 
'— ' , universality classes 0, Q ■ Recently, the average pulse shape has been proposed as a sharper tool for discriminating 
^ ' among universality classes and to test models against experiments 0. This analysis has revealed some weaknesses 
^ ' of present models, since they all fail to reproduce the avalanche shapes observed experimentally. Namely, all models 
ly-j proposed so far produce symmetric shapes, while leftward skewed forms are observed in experiments, indicating 
QQ that our understanding of the Barkhausen effect is, at the present stage, incomplete. This open issue has been the 
04 1 inspiration of this work. 

' We consider the problem of finding the average shape of a generic stochastic signal during an excursion, that is 
between two successive returns to a reference value. We hope that a deeper understanding of how the statistical 
properties of the signal are reflected on the shape of the average excursion can in general give insight into the 
understanding of the process generating the signal. In the case of Barkhausen noise, this may help identifying which 
"j^ • crucial ingredient is missing in the theory, and lead to the introduction of more accurate models. 

Beyond its interest for what concerns the understanding of Barkhausen noise, the nontrivial phenomenology of the 
I . avalanche shape leads to more general and interesting questions: What are the physical ingredients that determine the 
' shape of the average excursion in a generic stochastic process? Are there universality classes? Does this shape encode 
Ch pieces of information about the underlying physical system, that are not accessible by considering other observables? 
O These issues have not been addressed so far. In this paper we begin a systematic investigation of the shape of 
the average excursion, by considering some simple stochastic processes, both uncorrelated and correlated, and with 
generically distributed increments. In this way we provide a first theoretical framework that may be of help in the 
analysis of real time series in many contexts. 

In Section ^ we introduce the concept of excursion, the types of processes that we will consider in the following 
' and the general scaling form of the average excursion. Section UTTl presents the results for processes with uncorrelated 
increments (Brownian motion. Random walk. Levy flights) and Section llVI discusses the effect of a damping term in 
a Brownian motion. Sections IVl and IVTI consider, respectively, the effect of long- and short-ranged noise correlations. 
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FIG. 1; Schematic representation of the average shape of a fluctuation. Thin sohd hnes are two reahzations of the stochastic 
process x{t), both returning for the first time at zero at time T. The thick sohd hne is the average shape computed over many 
reahzations. 

The final Section IVIII presents some conclusions and an outlook. A short account of some of the results presented 
here appeared in Ref. 



Let us first define the average excursion of a stochastic process. We consider a real valued Id process x{t) defined by 
a Langevin equation with suitable initial conditions. An excursion of the process is the trajectory in the {t,x) plane, 
followed until the first return to the initial value a;(0) [see Fig. (^]. We are interested in the statistics of positive 
excursions of a given duration T, i. e. those such that x{t) > x{0) — for < t <T. In particular we will denote the 
average excursion as {x{t))T- 

When analyzing real experimental data, one may need to extend our definition to a generic reference value dif- 
ferent from zero. For example, this is the case for positive signals, for which the identification of excursion (avalanches) 
is made unclear by the presence of background noise: one has to choose a recipe to decide when an avalanche starts 
or ends. In practice, one sets a threshold which is small enough not to change the shape of the avalanche, but high 
enough with respect to background noise. 

A generic reference value a is taken into account by a translation of the origin of t and x: the excursion of x{t) with 
respect to a is the excursion of a process x'{t') — x{t' + ta) — a with respect to x'{0) = 0, where x{ta) — a. 

Notice that in general the probability distribution of first return times P{T) depends on the threshold value a. For 
this reason we expect that variations of a will affect the average excursion {x{t))T- In fact, for most of the simple 
processes discussed below the choice of a has no effect on the shape of the average fluctuation. 

We will analyze two kinds of processes. The first are processes of the type 



with random increments ^(t) extracted from a distribution Q{£,). In Section |IIII the noise ^ will be taken to be 
uncorrelated while in Sections IVl and IVTI we will deal with the effect of correlations. The case considered in Section Hvl 
is the simplest instance of a wide class of processes, Brownian motions in a potential 



For analytical calculations it is useful to express the average excursion more explicitly. Let us first introduce the 
excursion distribution Q{x,t \ xq, 0; xq, T), which is the probability that a trajectory, started in xq for t = and 
returning to xq for the first time at time T, is in x at time t. For each time t, is the distribution of the quantity 
whose average is {x{t))T, i. e. 



II. DEFINITION OF THE AVERAGE EXCURSION 



dtx{t) = m 



(1) 



dtx{t) 



dV{x)/dx + ^{t) 



(2) 



{x{t))T 



dx xQ{x,t I a;o,0;a;o,T) 



(3) 



dxQ{x,t I a;o,0;a;o,T) 
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Note that i}{x, t \ xq, 0; xq, T) is related to the distribution of first return times of the process P{T) = dx i}{x, t \ 

In the case of Markovian processes, as Eq. |2l, ^ can be written in terms of the function c{x, t \ xq, io)i which is the 
probabihty that the process started at xq at time is in x at time t, with the condition that x > for all to < t' < t. 
The probability fl(x,t \ xq,0; xo,T) of the whole trajectory is the product of the probability c{x,t \ a:o,0) of going 
from xq to x{t) times the probability c{xo,T \ x,t) to go from x{t) back to xq at time T. Notice that instead of 
starting exactly from we have to consider the path starting and arriving in xq ^ since c vanishes identically for 
Xq = 0. Hence, for any Markovian process, we can write 

J^dxc{x,t\xo,0)xc{xo,T\x,t) 

{x{t))T = Im -joo-; — ^-n ^71 7^ • (4) 

2:0^0+ ax c{x,t \ xo,U) c{xo, 1 \x,t) 

Eq. together with translational invariance and the scaling assumption 

c{x,t\xo,0)^t>^h{{x~xa)/t") (5) 

for the conditional probability c(x,t | xo,0) implies 

imr = TV{t/T). (6) 

In some of the cases that we will consider Eq. (@J cannot be applied, since noise correlations break the Markovian 
property. Nevertheless we will always find Eq. © to be true provided that the distribution P{T) of first return 
times decays algebraically. In all the cases considered the exponent a coincides with the wandering exponent of the 
unconstrained process, defined by {[x{t) — x(0)]^) ~ t". 

III. UNCORRELATED PROCESSES 
A. Brownian motion 

The simplest process is the uncorrelated Brownian motion 

dtx{t) = m (7) 

where ^{t) is a Gaussian white noise with {^{t)) — and {^{t)^{t')) = 6{t — t'). 

The probability c{x,t \ x^^Q) can be computed via the image method as a linear combination of two solutions 
of the Fokker-Planck equation associated to the free process 

C(X, t I Xo, 0) = ^ L-(x-xo)V(2t) _ g-(x+.o)V(2t)l ^ (8) 



yielding, in the limit of small xq 



The time-reversal and time-translational invariance of the process imply c{xq, T \ x^t) = c{x, T — t \ xq, 0), hence the 
distribution 17 is 

n{x, i I xo, 0; xo, T) = c(x, t \ xq, 0)c(x, T - i | xq, 0) cx [(T - t)t]-3/2(^2:o)'e-"'{i/(2t)+i/[2(T-t)]} (^q) 
Expression (O, inserted into Eq. Q), gives for the average excursion 



= W)T=W^A^fl-;). (11) 



n\] T \ T 

The average excursion of Brownian motion is thus of the scaling form JBJ, with the exponent a = 1/2 coinciding with 
the wandering exponent of the free process and a scaling function proportional to a semicircle 



fu{s)^\l-yMl^, (12) 
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where s = t/T. This result had already been noticed by Fisher Q- The variance of the excursion is also readily 
computed 

{ix-(xit))Tf)T^T(^3-iys{l-s). (13) 

The previous results are easily generalized to the case of a Brownian motion with bias, that is Eq. Q with 
V = > 0. The process is now invariant under time-reversal only provided the velocity is also reversed. Therefore 
in this case c(xo, T \ x,t;v) — c(x, T — t \ xq, 0; —v). Such quantities can again be computed via the image method Q 

c(x, t I xo, 0; v) = _l^e-("-"«-"*)'/('*) (l - e'"""/*) . (14) 

Inserting this expression into the formula for the excursion distribution it turns out that O is the same of the unbiased 
case except for an additional factor e"" Such a constant appears both in the numerator and the denominator of 
Eq. (0J implying that the average excursion is exactly the same of the unbiased Brownian motion 111|) . 

Notice that the addition of a bias introduces a characteristic time of order 1/w^, which reflects in a cutoff in the 
distribution of first return times ^ 

P{T) oc r-3/2e-'^/2. (15) 

However, a bias does not alter the shape of the excursion: the number of trajectories that survive up to a time 
T ^ is exponentially small, but the average shape of these unlikely events is exactly the same as for the 

unbiased case. 

This is the first example of the insensitivity of {x{t))T with respect to changes in the distribution of single steps, a 
feature that will turn out to be quite generic. 



B. Random walk 



The Brownian motion is a continuous process in space and time. On the basis of the central limit theorem it is 
reasonable to expect the form of the average excursion to be the same for all processes with finite variance of the 
single increments. To support this conjecture we now compute {x{t))T for a process with finite variance, discrete in 
space and time, a random walk with bimodal distribution of the noise, i. e. Q{^) = (^j,i + (5^,_i)/2. The number of 
paths starting in at time and ending in x at time t without ever touching the a; = axis is given by |^ 

Hence the probability to find the walker in x at time t with the condition that it has never touched the axis is obtained 
from Eq. 1)16(1 dividing F{x, t) by M(t) = -^(^i 0) the total number of possible trajectories of t steps in the positive 
X half-plane 

Using time reversal symmetry, the average excursion is then given by 

<^<'»- - ii"!:;^!:^,) = t " ^ - " ■ 

with K ^ M {t)M (T - t) /F{1, T - 1), where we used the fact that Y.x ^i^^ t)F{x, T - t) is independent from t and 
equal to F{1,T — 1). Introducing the variables 7 = aj/Vi, s = t/T, and = ^ s/{\ — s), and using the expansion 

n!/(n/2)! ^ y^2"/n, we get M{t)c{jVt,t) = F{jVt,t) ^7^-^'/^ i^(l,T- 1) ^ and K ^ 

so that the average excursion is 

{x{t))T = i^t^^^d77^e-*(i+^=) ^ (19) 
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FIG. 2: The factor N{T) vs T for Levy flights with several values of /i. Notice that the exponent is 1/2 for ^ > 2. 



where in the last step the integral has been extended to infinity. Expressing </> in terms of s ~ t/T we recover 



{x{t)), 



^^(^ 



(20) 



exactly as in the continuous case. 

In the case of a biased walk, where Q{S) — {qS^.i + (1 ^ one must replace the number of trajectories 

F{x, t) with their probability: each trajectory reaching x at time t is weig hted by a factor g(*+^)/2(i _ g)(t-^)/2. The 
trajectory leading back to zero has a instead a weight g(^~*~^)/2(l — q^(T-t+x)/2 ^ rpj_^g product of these weights gives 
simply a constant factor [q{l — q)]^^"^ both in the numerator and in the denominator of formula (|18() . Hence the 
introduction of a bias does not change all moments of the distribution of the excursion. 



C. Unbiased Levy flights 

Levy flights are statistical processes of the type ((Tjl, where the distribution of single steps Q(C) has a fat tail 
decaying as |CI ''^^ with < /j, < 2 so that their variance is infinite. The standard form of the central limit theorem 
does not hold for Levy flights: the invariant distributions under summation arc the Levy stable distributions flfil |. It 
is therefore natural to wonder whether Levy flights belong to a different universality class also with respect to the 
average excursion. 

The analytical evaluation of {xit))^ for this case is not straightforward as for the processes discussed so far. While 
Eq. HI) still holds, the image method cannot be used to determine c(x,t | a;o,0), because the Fokker-Planck equation 
for the free process is not local [T]| . We have therefore computed the average excursion numerically, considering steps 
performed at discrete integer times with absolute value distributed according to Q(|CI) « (|^| + 1)^^^^ and random 
sign. 

Details about the evaluation of {x{t))T in this and in the other cases where numerical results have been obtained, 
are presented in the Appendix A. 

To check the validity of the scaling hypothesis we compute the quantity N{T), such that {x{t))T = N{T)f{t/T), 
where we choose, with no loss of generality to normalize the scaling function / so that f{s)ds = 1. If scaling holds, 
N{T) has to be proportional to T" and the normalized average shapes {x{t))T /N{T) must collapse on the same curve 
for different T. In Fig. [21 we plot N{T) for several values of /x, showing that a — max[l/2, Also in this case a 

coincides with the wandering exponent of the free process, which is l//i for fi < 2 and is the usual diffusive one when 
^ > 2 12]. 

Fig. 13 reports the shape of the average excursion for values of /i such that the second or even the first moment of 
the single step distribution is infinite. In all cases the curves for different values of T collapse and the scaling form 
is exactly the same of the Brownian motion: The shape of the average excursion is completely independent from the 
distribution of single steps. 




0.2 0.4 0.6 0.8 1 



FIG. 3: Scaling function / for for Levy flights with several values of fi. 




10 



— t=100 

— t=1000 

IP 5P 

f(s)= a s /[ 1 + s ] 
a=5/(47i) sqrt(10)/sqrt(5-sqrt(5)) 



10' 



10 10 

s=x/t 



10 



10 



FIG. 4: Comparison between c{x, t) evaluated numerically for for Levy flights with steps distributed according to Eq. H21|l and 
the ansatz 12211 . 



In the Appendix B we report the results also for the variance of Levy flights. We have not been able to proof this 
result analytically for Levy flights. However, in the case = 1, we have considered the Levy-stable distribution, the 
Cauchy distribution 



Q{0 



1 



(21) 



We have computed numerically the probability density c{x, t) for such a process. The result is presented in Fig. 01 
where it is compared with the ansatz 



c{x, t) = — 



Heft 



(22) 



where a is a normalization constant. Formula (|22|) is the simplest expression that interpolates between the small s 
and large s power law behaviors found numerically. The agreement between the numerical results and the formula is 
striking. 

Inserting the expression l|22|l into Eq. Q and performing the integrals, one obtains 



{x{t))T = TVs(l-s). 



(23) 
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FIG. 5: Bias-dominated regime for Levy-distributed increments {fi = 1.5) with bias v — 20. Main: Normalized scaling function 
/ converging toward the asymptotic form f{s) — 2s. Upper inset: The factor N{T) growing as with a = 1. Lower inset: 
First return distribution P{T) decaying as T^**^^. 



adding another piece of evidence to the universahty of /(s). 



D. Biased Levy flights 

We now consider the case of biased Levy flights, where single increments are distributed symmetrically as in the 
unbiased case plus a constant term v. In this case the average form of the excursion is in general asymmetric for finite 
times (toward left or right depending on the sign of v). To understand this behavior, it is fundamental to consider the 
relative importance in the equation of motion of the drift term and of the wandering due to the stochastic variable ^. 

The former is clearly vT, while the latter grows as T^/^. Hence a crossover time T*{v) ^ y-t^/i'^-i-'-) exists between 
two regimes: which of the two mechanisms dominates depends on the value of /i. For 1 < /i < 2 the bias dominates 
for large times, while the wandering is larger than vT for T <tiT* . For /i < 1 the opposite is true and asymptotically 
the bias does not play any role. 

When the bias is irrelevant, the behavior is the same of the unbiased case: {x{t))T is given by Eq. H12|l . a = 1/^ 
and the first passage time distribution P{T) decays as y-i-i/A" [i^. 

When the bias dominates, instead, the shape of the excursion is completely different. In this case the trajectories 
are practically deterministic, i. e. ballistic motions with velocity v. However the noise term is crucial to have the 
constraint x{T) — satisfied, since at t ~ T a very large fluctuation is needed. As a consequence, the distribution 
of return times is P{T) ^ y-i-/^^ ^jjg exponent a is 1 and the average excursion has a triangular shape /t(s) = 2s. 
This bias-dominated regime is shown for /i = 1.5 and T ^ T* in Fig.|S| The crossover between the two asymptotic 
regimes is very slow and it is not possible to run a single simulation long enough to exhibit the full transition between 
the early and late regimes. Nevertheless it is possible to distinguish clearly between the behavior for /i > 1 and /i < 1. 
In the former case (Fig. 0) the form of the average excursion becomes more and more skewed with time, while in the 
latter case (Fig.[7|) the opposite behavior is observed. To give a measure of the asymmetry, we consider the quantity 

/ ds/(s) sign(s-l/2) (24) 
Jo 

which is zero for the semicircle and 1/2 for the triangle. The insets of Figs. El and [3 clearly indicate that, for fi = 1.5, 
it grows, while it decreases to zero for fi = 0.5, consistently with the argument presented above. 

Therefore, we can conclude that, while in the unbiased case the average excursion of a generic uncorrelated stochastic 
process with symmetric steps obeys asymptotically the scaling form © with universal shape /[/(s) [Eq. H12|) ]. in the 
biased case, this is true only for < 1 or /i > 2. For 1 < /i < 2 instead, the presence of a bias leads to an asymptotic 
average excursion with triangular shape. The results for biased flights are summarized in Fig. 151 
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FIG. 6: Main: Scaling function / for biased Levy flights with ^ — 1.5 and v = 1. Inset: Temporal evolution of the asymmetry 
parameter [Eq. (1241 1. 
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FIG. 7: Main: Scaling function / for biased Levy flights with fi = 0.5 and v — 25. Inset: Temporal evolution of the asymmetry 
parameter [Eq. 1241 1. 



IV. DAMPED BROWNIAN MOTION 



We now deal with a Brownian motion in a potential, Eq. (j^J. We treat only the simplest possible case, an harmonic 
potential V{x) = Ax^, so that the Brownian motion is pushed toward the origin by a linear damping term 

dtx{t) = -^Xx + e(t) (25) 

The formal solution of this equation is given by 

x{t) = xoe-^' + f dse"^(*-")?7(s) (26) 
Jo 

Since x{t) is linearly related to rjit), it has a Gaussian distribution characterized by its first and second moment, that 
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FIG. 8: Sketch of the average shape for positively biased flights. The crossover time T*, depending on the value of the bias, 
changes the asymptotic behavior only for 1 < fi < 2. For /i < 1 the average trajectory is asymmetric in the intermediate regime 
T <^T* , and it takes a triangular shape in the limit T oo after T* — > cxd. 



can be obtained by averaging over the noise: 



m{t) = {x{t)) = xo 



,-At 



The normalized probability density for the free process is therefore 



P(x,t) 



1 



■ exp 



(X - XqC'^*)'^ 



One can easily check that P{x,t) solves the Fokker-Planck equation 

1 



dtPix,t) = ^dlP{x,t) - \d^[xP{x,t)] 



associated to the Langevin equation H25|l. 
The distribution of first return times is 



1 p2AT 1 

P(r)cx— ^— - = — (2A)3/2 



-XT 



(1-e 



-2XTy/2 



Thus 



P(T) oc I i^"''' ^ « 

^^-^)^\ |(2A)3/2e-^i^ r»l/A, 



(27) 



(28) 



(29) 



(30) 



(31) 



For T <^ 1/A, the T"^/^ behavior characteristic of Brownian motion is recovered. When T is of order 1/A, the 
power law behavior is cut off exponentially. 

To calculate the average excursion we need to evaluate the probability c{x,t \ xo,0) that a walk originating at xq 
at time is found in x at time t, without having ever touched the origin. Since Equation 12511 is linear, the image 
method can be applied, yielding 



c{x,t I a;o,0) 



1 



2TTaxit) 



exp 



X — xoe 



-xt\ 



exp 



2a^x{t) 



In the limit xq —>■ 



-At 



/ , I n\ 2xQXe 2 It) 

c(x,t I xo,0) oc , „ , — e 



27ra3(t) 



(32) 



(33) 
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FIG. 9: Average excursion for a damped random walk [Eq. 13811 1 for 1/A = 20. From top to bottom lines are for T = 
10, 25, 50, 100, 250. Notice that the shape flattens to the constant value ^ 4/(7rA) (thin line). 



The quantity c(a;o, T | x, i) for — > is obtained in a similar way 

c(a;o, J I X, rj cx 

Defining ci\(t) = cr^(i)e^^* one can rewrite (j34|l as 



27ra3(r-t) 



OxnXe^^t^"*) '7^ 

c(xo, i I X, t) cx r- „ e a' ' 



(34) 



(35) 



Note that d\(t) = (T^(f)e^^* = cr^^(i), thus the process with reversed time formally corresponds to the process with 
-A. 

The distribution of excursions Q is then 



nix,t I xo,0;xo,r) cx e-^^[aA(t)a_A(r-0]-'e--'/P-^.(*'^)l 
where the variance of the Gaussian factor is 



1 

2A 



(l-e^2At)(i_g-2A(r-t)) 



1 - e 



-2AT 



(36) 



(37) 



Inserting Eqs. and into Eq. igji yields 



1 /(l-e-2>*)(l-e-2>('r-*)) 



TT ^/2A 



1 -e 



-2AT 



(38) 



As expected the existence of a characteristic time in the problem (1/A) breaks down scaling; the shape of the 
average excursion changes with T (Fig.^. However formula © still holds in the two asymptotic limits 



{x{t))T = 



fV^^^ i,r-t«i/A 

t,T-t > 1/A. 



1 



(39) 
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Thus in the small A limit the semicircle is recovered, while in the large A limit the curve flattens around a value 
proportional to 1/\/A, while keeping the semicircular tails. The crossover between the two regimes corresponds to a 
change in the variance (Jg^ of fl 

a' =r-^ t^T-t«l/X (40) 
I 21 i,T-t»l/A, ^ ^ 

Thus fluctuations saturate at the value 1/(2 A), indicating the absence of correlations on time scales longer than the 
characteristic time 1/A. 

In the computation presented here we have taken the reference value to coincide with the origin, i. e. a = 0. In the 
presence of a potential, the effect of a reference value different from zero corresponds to the change V{x) V{x + a). 
For the damped Brownian motion, this is equivalent to the introduction of a bias Xa. In general, such a bias perturbs 
the form of the average excursion, at odds with the case of free random walks, which are unchanged by the presence 
of a bias. However the change in the average shape is small provided a is small compared with A/^/ttX, the maximal 
amplitude of {x{t))T- 



V. LONG-RANGED CORRELATIONS 



We now start considering the effect of the introduction of temporal correlations in the process. In particular, 
we study processes of the form dtx{t) — £,{t), with correlations between single increments g{t,t') = {^it)£,{t')) — 

We first focus on a process dtx{t) — £^{t) where the noise performs in its turn a Brownian motion dt^{t) — rj{t) with 
(77) — and {rf) = 1. Clearly this process is non-Markovian and can be written as the Random Accelerated Particle 
(RAP) 

d^x{t)=r^it). (41) 

The correlation function of £,{t), g(t,t') = mm{t,t') does not decay to zero when t — t' diverges: the noise has then 
infinitely-ranged correlations. This process has been studied recently with regards to polymers [Tsl l and the inelastic 
collapse of granular matter . 

It is important to stress here that, at odds with the previous cases, the non-Markovian nature of the process implies 
that to define completely a fluctuation one has also to consider the initial and flnal velocities (the process is Markovian 
if one considers the broader space of coordinate and velocity). We consider separately avalanches beginning with zero 
or a finite velocity, that will yield different results. For what concerns the final velocity the most natural thing is to 
average over all final velocities. The condition on the initial velocity is numerically very easy to implement, since it 
corresponds to setting 77(0) = vq. Instead, the condition on the final velocity is more delicate. Considering returns 
within a strip (— e, +e) implies, when using discrete times, a hidden condition on the final velocities. This is apparent 
if one considers the long time decay of the distribution of first return times. If one considers discrete times and a 
return in a strip, the known result P{T) ~ 2.7] is not recovered. The expected exponent is found instead if one 

averages over all trajectories positive up to t < T and becoming negative for t = T . Therefore we considered only this 
latter kind of trajectories. We believe that, asymptotically, this coincides with the continuous time process averaging 
over all final velocities. 

Figure [TUI reports the average excursion shape for two different distributions of rj with finite variance, one uniform 
and the other exponential, both in the case of zero initial velocity. It is clear that scaling holds very well for all values 
of T considered and for both distributions of the noise, and that the scaling function is asymmetric. As shown in the 
inset, the exponent a is again equal to the free wandering exponent, which is 3/2. Remarkably, the scaling function 
is indistinguishable from the simple form 

f{s)=^-s^'\l^s). (42) 

In the Appendix B we report also the results for the variance of the excursion, showing that it scales as the first 
moment. 

When the initial velocity is finite, instead, the asymptotic scaling associated with an average shape H42II is 
preceded by a transient regime with a more symmetric {x{t))x, depending on T and Vq (Fig. [TT)l. For large vq the 
shape is very close to a simple parabola 6s (1 — s). We have also checked that the exponent a crosses over from a 
ballistic value a = 1 at short times to the asymptotic value 3/2. Comparing the values of a in the two regimes we 
expect the crossover time to be proportional to Vq. 
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— Eq. 42 
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FIG. 10: Main: Average excursion for a RAP with zero initial velocity, with uniform or exponential distribution of the noise ri, 
showing perfect agreement with the simple form 14211 . Inset: Factor N{T) and first-return time distribution P{T) for uniform 
noise distribution. 




FIG. 11: Main: Scaling function / for a RAP with finite initial velocity, with uniform distribution of the noise rj. 

We then consider RAP with noise rj distributed with slowly decaying tails i^(|?7|) ~ rj^^^^. With zero initial 
velocity the scaling function depends on the value of (Figure I12|l . at odds with what occurs for uncorrelated 
unbiased processes. The skewness is toward right, the more so for small values fi. The exponent a is equal to 1 + 1/^. 

Despite the dependence of the detailed scaling form of {x{t))T on the distribution of single steps, some degree of 
universality remains also for the RAP process. The right panel of Figure El shows that the exponent characterizing 
the behavior of /(s) for s 1 is universal, being linear also for fi = 1. For small values of s instead, a careful analysis 
indicates that the actual exponent of the power-law behavior is slightly different from the value 3/2 obtained for finite 
variance (Figure ^1 left panel). We do not have an explanation for the value of this exponent. However, we cannot 
rule out that such universality is restored for larger values of T. We have checked that, in the case of finite initial 
velocity vq^ the transient parabolic shapes are present also for values of /x < 2. 
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Eq. 42 




FIG. 12: Main: Scaling function / for for a RAP with noise rj distributed with slow decaying tails P{\ rj |) and several values 
of fi, compared with the simple form 14211 . Inset: Factor N{T) and first-return time distribution P(T) for the same values of /i. 
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FIG. 13: Left top: Small t/T tail of the scaling function / for RAP with Cauchy distribution of single steps. Left bottom: 
Local effective exponent computed on the figure above. Right top: Small 1 — t/T tail of the scaling function / for RAP with 
Cauchy distribution of single steps. Right bottom: Local effective exponent computed on the figure above. 



VI. SHORT-RANGED CORRELATIONS 



We now turn to the case of short-ranged correlations. We consider a process with correlations decaying exponentially 
over an interval r 

g{t,t') = e^p{-\t-t'\/T). (43) 

A process of this type can be obtained in practice by feeding the random walk with noise obeying a damped random 
walk of the form + 1) = j£,it) + — j'^)'ri{t) with uncorrelated rj. The correlator of ^(i) is easily shown to decay 
exponentially, with a characteristic time r = — l/log7. 

For such a process, we do not expect scaling to hold for all values of T. We can anticipate instead two regimes 
depending on the duration T of the trajectories considered. For short times T <ti t, noise is correlated during the 
whole trajectory and the behavior must be the same of the case with infinitely-ranged correlations treated above. For 
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FIG. 14: Normalized average shape for a short-memory process with zero initial velocity, finite variance noise and t = 1000. 
The dashed lines are the expected limiting curves for T <C r and T ^ r. 

long times T ^ r, on the contrary, noise is correlated only for intervals that are short compared to the total duration 
of the excursion. Hence the process is equivalent to an uncorrelated process with some effective distribution of the 
single increments. 

Numerical results fully confirm this picture. In Fig. 1141 vye show the case of a short-memory process with zero initial 
velocity and finite variance noise. For short times T <C t the shape is very close to the form (|42() . valid for RAP. A 
slow crossover leads for longer times to the semicircle law valid for uncorrelated processes. We have checked that the 
expected pattern of behavior occurs also for non zero initial velocity and for Levy-distributed noise. 

VII. CONCLUSIONS 

Let us summarize the results presented in the previous sections. We have studied the statistics of excursions in 
some classes of stochastic processes, with particular attention to the average shape. For uncorrelated free processes 
we have found that the average excursion has a scaling function proportional to a semicircle, independent from the 
distribution of the single steps provided it is symmetric. This holds not only for distributions that renormalize to 
the Gaussian, but also for the class of distributions that renormalize to symmetric Levy stable distributions. More 
generally the scaling function is unchanged when a bias is introduced, with the notable exception of the case with 
1 < fj, < 2, where the asymptotic shape of fluctuations is triangular. The addition of a linear damping term in the 
Langevin equation for the process introduces a characteristic time scale, that separates between two regimes: for short 
times the process is dominated by noise, and the excursion is the same as in the free case. For longer times, scaling 
breaks down and the shape of the average fluctuation flattens to a value independent from its duration. Furthermore, 
we have analyzed the effect of noise correlations for the free process. When correlations are long ranged the shape of 
a fluctuation depends on the initial velocity ■^o- For vq = we find that the scaling function has asymmetric tails s'^/^ 
and (1 — s) in the Gaussian case, while the situation for Levy distributed steps is less clear. For vq > a transient 
regime exists such that the scaling function has linear tails (independent from the distribution of the single steps), 
before it crosses over to the asymptotic form which is the same of the vq — case. Finally, in the case of short 
range correlated noise, the range of correlation sets a time scale that separates between a short time regime, where, 
as expected the behavior is similar to the long range case, followed by a crossover to the asymptotic uncorrelated 
behavior. 

Application of this analysis to real data requires some care. Indeed in many situations of practical interest one 
deals with long time series consisting of a large number of successive fluctuations. In such a case, if one computes 
{x[t))T by averaging over successive returns to the value x = a one may average over pulses that are not statistically 
independent. 

In our work, on the contrary, we take care to average always on independent events. When the process we consider is 
Markovian this does not require particular prescriptions. In this case averaging over successive fluctuations in a single 
realization is equivalent to averaging over avalanches belonging to different realizations. Otherwise one should consider 
avalanches separated by times larger than the largest correlation time in the system, if the analysis is restricted to a 
single realization. If the correlation time is inflnite, as in the RAP, one should consider only fluctuations belonging 
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to independent realizations. 

Another relevant issue for the application to real time series concerns the amount of events required to obtain 
sufficiently clean results. In principle one should average over fluctuations of exactly the same duration, and rescale 
afterward. This may turn out to require an exceedingly large number of events. An alternative procedure is to assume 
scaling and average over fluctuations of different duration, properly rescaled, with an exponent that can be obtained 
by plotting the size of fluctuations as a function of their duration. This also checks whether scaling holds or not. 

For the case of Barkhausen noise, which was the initial inspiration of this work, we suspect that the asymmetric 
shape observed in experiments must be due to the presence of some kind of correlations. However, the kind of 
correlations that we have analyzed give rightward asymmetric shapes, while the one observed experimentally are 
leftward. This calls for further analysis, of more general processes. 

APPENDIX A: DETAILS ABOUT THE NUMERICAL RESULTS 

When performing numerical simulations, we have taken times to be discrete and space continuous, thus the concept 
of first return to the initial value needs some clarification. We have considered the process to return to the initial 
value when its value is in a small interval around it [— e, e]. In order not to introduce an artificial asymmetry we 
have applied the same condition to the first step of the excursion as well. This is implemented by letting the process 
start at a; = for some negative time and taking as i = 1 the first time such that x(t) > e. The average is then 
performed over all trajectories that first return between — e and e at a specified time T under the constraint that 
x{t) > e for 1 < t < r. Care has to be used when choosing the value of e, which should be as small as possible to give 
e independent results. At the same time too small values of e make the numerical simulation very time consuming. 

In all simulations T is integer; hence it is in principle possible to average only over trajectories that return after 
exactly T steps. However T cannot be decided a priori but is the outcome of the simulation; collecting a large number 
of trajectories for a single large T may be a prohibitive task. Therefore, similarly to what is done in experiments, we 
have performed a binning procedure, by averaging over all runs that return within a time interval between T — AT 
and T + AT. Since the relative size of the bins 2AT/r decreases as T grows, we are sure that this procedure does 
not lead to significant perturbation of the results for large T. 

Finally, when presenting numerical evaluations of the average excursion we plot {x{t))T normalized by the factor 
N{T) = Jg ds{x{sT))'j'- This allows to check the form © in two ways. If scaling holds then N{T) grows as a 
power-law (the exponent is a) and the different curves collapse on a universal shape [which is /(s)]. 

APPENDIX B: NUMERICAL COMPUTATION OF VARIANCE 

For the simplest processes that can be attacked analytically we have already computed the variance: 

4(t) = {ix{t) - {x{tm 

In the other cases it is possible to determine numerically such a quantity. In this appendix we show, as an example, 
the results for two of such cases, i.e. the unbiased Levy fiight with fi — 1.5 and the RAP. In Figures [T^ and [TBI 
we have reported the average excursion {x{t))T for different durations T divided by the expected scaling factor T". 
Correspondingly we plot the values of the standard deviation o'T(i) divided by the same factor. One sees that a scales 
exactly as the average shape in both cases. 
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FIG. 15: Plot of the average shape and standard deviation for Levy flight with jj. — 1.5. Curves for different durations are 
divided by the scahng factor T^^>^ . The upper curves are the average shapes, the lower ones are the standard deviations. 
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FIG. 16: Plot of the average shape and standard deviation for RAP with Gaussian noise. Curves for different durations are 
divided by the scaling factor T^/^. The upper curves are the average shapes, the lower ones are the standard deviations. 
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